Oscillatory spin relaxation rates in quantum dots 
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Phonon- induced spin relaxation rates in quantum dots are studied as function of in-plane and 
perpendicular magnetic fields, temperature and electric field, for different dot sizes. We consider 
Rashba and Dresselhaus spin-orbit mixing in different materials, and show how Zeeman sublevels 
can relax via piezoelectric and deformation potential coupling to acoustic phonons. We find that 
strong lateral and vertical confinements may induce minima in the rates at particular values of the 
magnetic field, where spin relaxation times can reach even seconds. We obtain good agreement with 
experimental findings in GaAs quantum dots. 
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Since the proposal of a qubit based on the electron spin 
of quantum dots (QDs),i much work has been done to un- 
derstand the processes that may cause their relaxation, 
since long coherence times are required. One of those 
processes is related to the phonon-induced spin-flip rates 
of Zeeman sublevels in QDs in magnetic fields, where the 
spin purity of the levels is broken by the spin-orbit (SO) 
interaction. A recent experiment^ has shown that spin 
relaxation time has a lower bound of 50 fis at an in-plane 
field of 7.5 T in a GaAs QD defined in a 2DEG, while 
theoretical work'^ has given 5 /is at the same QD parame- 
ters. In general, SO effects have been considered via per- 
turbation theory,^"'' although exact treatments have also 
been presentedi^i^ The perturbative approach, which 
includes only a few states, has been called into question 
by the demonstration that a larger basis is needed in or- 
der to achieve convergence even for the lowest QD states 
when the QD vertical width is narrow^ where a complex 
interplay between different energy scales is presenti^ The 
influence of high SO coupling on the QD effective ij-factor 
has been recently addressed.^" 

Using exact diagonalization, our goal is to provide a 
bridge between the different approaches used so far, so 
that their range of validity can be well established in 
terms of materials (wide/narrow-gap), QD sizes, and di- 
rections of the magnetic field. For GaAs we show that the 
perturbative approach is indeed satisfactory for in-plane 
fields for any QD confinement, where good agreement is 
found with calculation2i& and experiment;^ and where the 
piezoelectric potential is shown to dominate. If the field 
is perpendicular, however, we find that the perturbative 
approach** indeed fails, giving relaxation rates orders of 
magnitude longer than our results if the z- confinement 
is strong^ For higher perpendicular fields, we show that 
minima are found in the rates depending on the QD fea- 
tures. For InSb the oscillatory rates are found for both 
in-plane and perpendicular fields and strong confinement. 
We find that the deformation potential dominates the 
rates, so that previous results^ under in-planc fields not 
including such potential yielded rates orders of magni- 
tude smaller than ours. We show that the oscillatory 



behavior is caused by the nature of the z-confinement, 
as well as by the complex interplay between the physi- 
cal scales playing a role in the system: Zeeman and SO 
energies, as well as QD lateral and vertical sizes. 

The QD is defined by an in-plane parabolic confine- 
ment, V{p) — mujQp'^/2, where m {ujq = Eq/K) is 
the electronic effective mass (confinement frequency); 
the QD lateral length is /q — x/fi/imujo)- The ver- 
tical confinement V{z) is strong enough so that only 
the state in the first quantum well subband is relevant, 
and its function is (fiziz) = y^2/zo sin (ttz/zq) if a hard 
wall is assumed, zq being the QD vertical well thick- 
ness. In a magnetic field B, the unperturbed Hamil- 
tonian, Hq = h^k^ /2m + V{p) + Hz, has the well-known 
Fock-Darwin (FD) solution,*-'-- where Hz — gofJ-B^ • c/2 
is the Zeeman term. We include all SO terms in 2D 
zincblende QDs, namely RashbaiS and Dresselhausi^ in- 
teractions. The former is due to surface inversion asym- 
metry (SIA) induced by the 2D confinement, while the 
latter is caused by the bulk inversion asymmetry (BIA) 
present in zincblende structures. The SIA Hamiltonian 
is HsiA = eta ■ Vy(p, z) x k, while the BIA is Hbia — 
iWrK m - kl)+ayky {kl - kl)+azkz {k^ - k^)], with 
coupling constants a and 7. The z-confinement yields 
the electric field dV/dz in the SIA Hamiltonian as well as 
the momentum average (k^) ~ ('"'/•^o)^ in the BIA terms. 
The full QD Hamiltonian is then H = Hq+Hsia+Hbia, 
which is diagonalized in a basis containing 110 FD states. 
Details about the derivation of terms in H, as well as 
their selection rules for (anti)crossings in the QD energy 
spectrum, are found elsewhere&ifi 

We calculate spin relaxation rates between the two 
lowest QD Zeeman sublevels caused by piezoelectric 
and deformation acoustic phonons via Fermi's Golden 



|'|M,(Q)|2(nQ- 



Rule: Tf, = 27r/?iE, q \Z{q. 
1)6 {AE + hcjQ), where the sum is over the emitted 
phonon modes j {j — LA,TA\,TA2) with momen- 
tum Q = (q,g,). The term Z(q,) = e^«-^ |^,) 
(7/i(ct) — (/I e**'''" N)) is the form factor perpendicular 
(parallel) to the 2D-plane (position is R = (r, z)), while 
nq is the phonon distribution with energy hcjQ; energies 
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AE = Sf—Ei and states \f) are obtained via diagonal- 
ization of the total H, so that the SO mixing is fully taken 
into account. The element M, (Q) = Aj(Q) +iSj(Q) in- 
cludes both piezoelectric Aj and deformation Sj poten- 
tials; in zincblende structures^ and in cylindrical coor- 
dinates, they become Ela{Q) — 'EqAlaV^ (only LA is 
present for S^), Ala{Q) = SAq^la sin(26»)g2gj2QV2^ 
Atai(Q) = AoATACOs(20)qqJQ5/2, and Ata2(Q) = 
AoATASm{29){2ql/q^ ~l)qy2Qy^ (both TAl and 
ryl2 modes are compacted as a single TA mode for 
Aj), where Aj = ^/h(2NoVc/y-^ and Aq = inehu/K. 
The bulk phonon constants are Sq and ehu, Cj are 
the sound velocities {ctai = CTyi2 = cta 7^ cla), k is 
the dielectric constant, and A^'o is the electron density. 
The triple space integration ([r, yields an ana- 

lytical solutionji and we were able to do two {[4>q,qz]) 
out of the momentum integration, leaving a numer- 
ical integral only in q. For a later use, the only 
zp-dependence in this remaining integral in Tfi reads 
Fj{zo) = {djZo — {djZo)^ /tt^)"^ sm'^{djZo), where dj = 
y/{AE/hcj)'^ - ^2/2; q runs from to AE/hcj, while 
Fj{za) is multiplied by polynomials and exponentials in 
q in the total T fi. No approximation is needed in our 
derivation of Fj^i, so that the 3D nature of the phonon is 
taken into account and the full form factor e'*^-^ is used. 

Panel A of Fig. ^ shows the phonon- induced spin-flip 
rates as a function of in-plane field, B\\ , for different val- 
ues of Eq (solid lines) for GaAs QDsii^ Larger i?o (smaller 
QDs) present smaller spin-flip rates, i.e., longer relax- 
ation times. This happens because as E^ increases, the 
orbital levels become more separated and then the SO 
coupling becomes relatively less important. A strong de- 
pendence of the rates with J5|| and Eq is clear. For ex- 
ample, at B\\ = 2 T, one gets 2.5 x 10^ s'^ at Eq = 0.7 
mcV and 2.5 x 10"^ s^^ at E^ = 5.0 mcV; at B\\ = 15 T 
the rate goes from 5.2 x 10^ s~^ to 6.9 x 10^ s~^ for those 
same values of Eq. Finite temperature (dashed-dotted 
line for T = 12 K) enhances the rates by about one or- 
der of magnitude due to the enhancement of the phonon 
population. Larger wells (dotted line for zq = 100 A) 
have rates about one order of magnitude smaller because 
the BIA term ((fc^) ~ l/^o) decreases. The same T and 
zo behavior is found for any Eq. Larger electric fields 
(dashed line for dV/dz = —38 meV/A) have little effect 
on the rates for large Eq, but at smaller Eq the rates 
decrease by one order of magnitude; this means that by 
increasing the SIA term so that it equals the BIA cou- 
pling, the rates can decrease even though the spin mixing 
gets stronger. For GaAs QDs under in-plane fields, we 
find that the TA piezoelectric coupling dominates at low 
fields (< 14 T for Eq = 5 meV), while at high fields the 
deformation potential takes over (symbols in Fig.^A). 

Panels B and C of Fig. ^ respectively show the spin 
expectation value (ct^) of the Zeeman sublevels and their 
energy splitting AE. Smaller Eq values show higher SO- 
induced spin mixing, so that (cr^) deviates from the pure 
-1-1 (—1) value for the ground (first excited) state, and 



FIG. 1: (A). Zeeman sublevel spin-flip rates for GaAs QDs— 
under in-plane magnetic field for different QD confinements 
(solid lines). Dashed-dotted, dotted, and dashed lines respec- 
tively show influence of higher temperature, well width, and 
electric field on the 3.5 meV QD. Symbols show isolated con- 
tributions from the distinct phonon mechanisms for the 5.0 
meV QD; TA piezoelectric dominates at low fields. (B). Spin 
average values of the Zeeman sublevels after diagonalization of 
H. Squares (circles) refer to ground (first excited) state. (C). 
Energy splitting of Zeeman sublevels after diagonalization of 
H. Dashed-dotted lines show results without SO coupling. 
Color scheme in panels B and C follow from panel A. 



FIG. 2: Same as Fig.0 but for GaAs QDs under perpendic- 
ular magnetic field. Vertical lines in panel B show the field 
where a sublevel crossing occurs and the normal spin charac- 
ter of the states is restored. Arrows in panel C indicate where 
suppression of SO coupling is induced by the magnetic field, 
which is at the origin of the respective minima in panel A. 



AE deviates from the pure qq^bB value of a QD with- 
out SO coupling. Notice at the 3.5 meV QD how a larger 
well (stronger field) decreases (increases) the SO-induced 
spin mixing by reducing (enhancing) the importance of 
the BIA (SIA) coupling. For comparison, the experimen- 
tal lower bound for the spin-flip time is 50 /iS at B\\ = 7.5 
T for a 1.1 meV QD,^ while from panel A we find 5 /is. 
As discussed above, this time is increased by one order of 
magnitude in a larger z-well, which appears to be a bet- 
ter match for the experimental conditions. A 5 /iS value 
has also been found from a perturbative formulation;^ 
and our results also agree with a different model for the 
vertical confinement^^ One can then affirm that the per- 
turbative approach^ is indeed enough for GaAs in in- 
plane fields at any QD sizes, and that the form of the 
z-confinement has little effect on the rates. 

In Fig. 121 we do the same analysis for GaAs QDs^^ in 
a perpendicular field, B^. The most striking feature is 
the appearance of minima in the rates shown in panel A. 
The £'o-dependent minima (12 < B± < 16 T) are due to 
the vanishing AE, as shown in panel C. At such values 
of B±, a level crossing of Zeeman sublevels results in a 
sudden spin-flip, as verified by the vertical lines of panel 
B for the 0.7 and 3.5 meV QDs. To understand the min- 
ima at low fields (6 < B± < 8 T) we have to be mindful 
of the sine argument in Fj(zo), which may induce a min- 
imum in the rate at a particular value of the field, due to 
the interplay of energy and length scales in the problem^ 
This behavior does not produce minima for in-plane field 
on GaAs QDs (Fig. . We find that those three low-field 
minima in Fig.[2J4 occur at B± values where a magnetic 
field-induced cancellation of the SO influence is produced 
on the respective AE values, as indicated by the arrows 
in panel C: below (above) such i5o-dependent values of 
B± - where AE recovers its Zeeman value of gofisB - 
the SO coupling increases (decreases) the Zeeman split- 
ting as compared to the QD without SO. 
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FIG. 3: InSb QDs^-*^ under in-plane magnetic field. The con- 
finements are such that the same QD size-range used for GaAs 
is covered. In panel A, the 15.0 (20.0) meV QD is chosen for 
the study of different QD parameters (distinct phonon mech- 
anisms); dashed arrows show that both LA modes have the 
same number of rate minima and happen at the same fields. 

Notice that the two smallest confinements do not show 
the low-field rate minima in Fig. |2H- This is due to the 
lowest sublevels acquiring the same spin (between 3 and 
7 T at i?o = 0.7 meV in panel B), so that the rate de- 
creases monotonically at that field-range. These features 
have strong influence on QD effective g-factorsiiS Like in 
the in-plane field problem, higher temperatures (T = 12 
K) enhance the rates by one order of magnitude, and 
the dominating phonon mechanism (symbols) is the TA 
piezoelectric coupling, but now at any field. To confirm 
the intricate interplay between all QD energy scales, a 
larger well (zq — 100 A) or a stronger field (dV/dz — —38 
meV/A) removes the rate minima, as shown in panel A 
for the 3.5 meV QD. Our results agree with available 
calculations^ at small fields {B± < 1 T) , so that we can 
confirm that the perturbative approach is not adequate 
when dealing with SO effects in QDs in a perpendicular 
field, even if the material is GaAs, if the vertical confine- 
ment is strong enoughi^ From Figs. ^ and |21 the rates at 
small fields (< 2 T) have the same behavior under 
and B±, namely, the smaller Eq the larger the rate; how- 
ever, they are much larger under B±. The rate maximum 
(at 2.5 T for the 1.1 meV QD) is due to an anticrossing 
involving the first and second QD excited statesiiS 

FigureOshows results for InSb QDs^^ in in-plane fields. 
In contrast to GaAs (Fig. , minima in the rates are vis- 
ible under B\^ for the highest Eq values. At small fields 
(< 2 T), a monotonic rate drop is again observed as Eq 
increases. As shown for the 20 meV QD, the LA defor- 
mation potential is the dominant phonon mechanism for 
InSb at any field (and Eq); this seems to be the case 
for any narrow-gap material. It is worth mentioning 
that, contrary to GaAs in Fig. [2H where a unique low- 
field minimum is present for a given Eq, the number of 
rate minima in InSb increases with Eq. Also, notice that 
both LA modes have the same number of minima at the 
same values (dashed arrows in panel A), while a larger 
number of minima in the TA mode occurs at different 
fields; this happens since Cla > cta, so that the sine 
argument in Fj{zo) is smaller in the LA mode. Panel B 
shows that Zeeman sublevels have normal spin character, 
where ground (excited) state is predominantly spin-up 
(spin-down) at low fields. Finite temperature (T = 12 
K) has no appreciable effect on the rates, a reflection of 
the much larger energy scales here. The energy scale in- 
terplay results at higher field {dV/dz = —2 meV/A) in 
slightly smaller (larger) AE (rates), and rate minimum 



shifts to higher fields. On the other hand, a larger well 
{zq — 100 A) produces a higher AE, so that the spin-fiip 
rate is suppressed by orders of magnitude and its num- 

FIG. 4: Same as Fig.|^ but for InSb QDs under perpendicular 
magnetic field. The only minimum of the 3.0 meV QD around 
1 T in panel A indicates a sublevel crossing, as shown by the 
vertical line in panel B and the vanishing of AE in panel C. 
The color scheme is the same in all panels and figures. 



ber of minima increases; notice in panel B that the two 
lowest levels acquire the same spin at ~ 8 T for such Zq, 
so that the oscillatory rate is well defined until that field. 
Previous calculations^ in InSb did not include the defor- 
mation potential, and yielded unrealistically low rates. 

For InSb QDs in perpendicular fields, we see from panel 
B in Fig. ^ that all spin-flip rates in panel A are mean- 
ingful only for B± < 8 T, a field where both Zeeman 
sublevels acquire the same spin. This happens earlier for 
small Eq as shown by the 3.0 meV QD, whose rates are 
well defined only for B± < 2 T; the vertical line around 
1 T in panel B for this Eq indicates a spin-flip, accompa- 
nied by the vanishing AE in panel C, the rate minimum 
in panel A, and a sign change of the QD factor Larger 
values of Eq present normal spin behavior. Spin-flip rates 
in panel A do not show the monotonic behavior seen for 
small fields in Fig.O Like in the case, a finite temper- 
ature (T = 12 K) does not show visible influence on the 
rates under B±, and the deformation potential dominates 
at any field. A higher electric field {dV/dz = —2 meV/A) 
slightly increases the rate; a larger well (zq — 100 A) in- 
troduces an oscillatory rate orders of magnitude smaller. 
We emphasize that the perturbative approach finds no 
use in InSb QDs because of the inherent higher SO cou- 
pling. It is worth mention that oscillatory rates have also 
been found in GaAs QDs under B± by considering mo- 
mentum relaxation from an excited orbital level to the 
ground state, as well as in coupled GaAs QDsi^ 

Even though both QD materials show oscillatory spin- 
flip rates - with B± for GaAs and with both _B|| and 
B± for InSb - their origin is slightly different. Minima 
come from the nature of the ^-confinement, and the field 
where they occur depend on the lateral size Iq. The SO 
coupling mixes spins and alters splitting of sublevels in 
distinct ways according to field-direction and QD mate- 
rial, so that spin relaxation can be induced by piezoelec- 
tric (wide-gap) and deformation (narrow-gap) phonons. 
Our calculations reveal the rich interplay between all rel- 
evant energy scales in QDs. The external field opens 
channels (at the rate minima) where long spin relaxation 
times (~ 1 s) may be reached, so that the spin coherence 
required for quantum computing could be improved. 

We thank support from NSF-IMC grant 0336431, 
CMSS at OU, and the 2P* Century Indiana Fund. 
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